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ABSTRACT 

Angular momen t um t ransfer from a rotating stellar bar has been proposed by 
IWein ber~ &Katz'("2002^ as a mechanism to destroy dark-matter cusps in a few rotation peri- 
ods. The A'-body simulations performed by these authors in support of their claim employed 
spherical harmonics for the force computation and were, as shown by Sellwood (2003), very 
sensitive to inclusion of asymmetric terms (odd /, m). In order to disentangle possible numer- 
ical artifacts due to the usage of spherical harmonics from genuine stellar dynamical effects, 
we performed similar experiments using a tree code and find that significant cusp destruction 
requires substantially more angular momentum than is realistically available. However, we 
find that the simplified model (a A^-body halo torqued by a rotating bar pinned to the origin) 
undergoes an instability in which the cusp moves away from the origin. In presence of this 
off'-centring, spherical density profiles centred on the origin dis play an apparent cusp-removal. 
We strongly suspect that it is this effect which Weinberg & Katz observed. When suppressing 
the artificial instability, cusp removal is very slow and requires much more angular momentum 
to be transferred to the halo than a realistic stellar bar possibly possesses. 

Key words: Methods: A^-body simulations - Galaxies: kinematics and dynamics - Galaxies: 
evolution - dark matter 



1 INTRODUCTION 

Cold Dark Matter (CDM) cosmological simulations consistently 
predict that the density profile of dark matter halos increases 
sharply towards the centre, creati ng a 'cusp', where the den sity in- 
creases like a inverse power-law. iNavarro. Frenk & Whitd il99dh 
proposed 



p(r): 



(r/r,)(l + r/r,)2 



(1) 



as a universal density profile (hereafter called 'NFW') for CDM 
haloes of all masses. Here, is a scale radius and is a characteris- 
tic density. This density profile has a central slope of d Inp/d In r = 
- 1 . Later work has questioned this inner slope and the universality 
over all halo masses, but the existence of a cusp with dlnp/dlnr 
in th e range -1 to -1.5 is a ro bust prediction of such simulations 
(e.g. lMoore et alll998l:rPower et aL2003 V 

Observations of the rotation curves of disk galaxies are the 
best probes of dark matter on galactic scales, and while some 
stud ies cannot rule o ut co nsistency with a cuspy NFW -like pro- 
file <Swaters et a l.'2003: van den Bosch & Swaterj200lh . the ma- 
jority of the observatio ns are claimed to be inc onsistent w ith the 
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'2001a"b"; 'Blais-Ouellette et al.' '2001: iSalucci et alJ 
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2004). Investigations combining rotation 



faavashi et alj2004 . it is considered the greatest current challenge 
of the CDM paradigm. 

Numerous possibilities have been put forward in an effort 
to explain the difference between theory and observation. Many 
of these involve altering the assumed properties of dark matter, 
either the dark matter particle itself I Soergel & Steinhardt 200ol 
iKaplinghatetalJbOOd) . or t he initial power spectrum on smal l 
scales (warm dark matter; e.g. lColfn et all2000l:[Bode et all200lh . 
Since the central region of most galaxies are dominated by bary- 
onic matter, it seems reasonable that the interaction between bary- 
onic and dark matter will have a major effect upon the dark matter 
density in the inner region of the halo. If the efi"ect of this interac- 
tion was to reduce the dark matter density at the centre of the halo, 
then this could potentially solve the problem. Several investigation 
of such processes have been undertaken, some pointing out effects 
that would increase the gradient of the density cusp still further 
I Blumenth al et alll986h . s ome discussing processes that could po- 
tentially remove the cusp fBinnev et aT]|200ltlELZant et alj|200lt 
[Weinberg & Katz 2002). In this paper we investigate the mecha- 
nism proposed bv (Weinberg & Katzi of angular momentum trans- 
fer from a rotating disc bar to the halo as a means of removing a 
cusp. 

JWeinbe^^^Cat^ po inted to an earlier study by 



poll 

iHernguist & Weinberg! \\99 2) to suggest that angular momentum 
would be rapidly transfered from a rotating disc bar t o the halo 
at an inner Lindblad-like resonance. IWeinberg & KatJ supported 
the hypothesis with analytical calculations and simplified A'-body 
simulations. In the simulations, the gravitational potential of a 



curves with gravitational lensing o bservations came to the same 
conclusion <Trott & Webste3l2002l) . While it has been suggested 
that the discrepancy may be a result of direct comparison with the 
NFW fitting formula, rather than with more realistic triaxial halos 
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rotating bar, with centre pinned at the initial centre of the density 
distribution, was imposed upon an A'-body representation of a 
CDM halo with an NFW profile. 

The simulations of IWeinberg & Katzl were re-examined by 
ISellwooj j2003) who found the same reduction in central den- 
sity, but argued that the deliberately artificial perturbation applied 
to the system lead to misleading results. Sellwood (2003) ran fully 
self-consistent simulations, which showed a slight steepening of 
the halo's inner profile from the action of a bar. This result was 
also seen in the simulations of Valenzuela & Klvpin ( 20 oJ). How- 
ever t he self consistent simulations of Hollev-Bockelmann et alj 
showed a c lear flattening in the cusp, as do the later ones 
of lWeinbergH2004) . A di scussion of this appa rently irreconcilable 
discrepancy was given by'Athanassoula' j 2004h . 

T he simulations of Weinbere & K at3 J20p2h and ISellwo oj 
jioO^), and also the original study of lHernauist & Weinberg ( 1992) 
all employed spherical harmonics to compute the gravitational po- 
tential and accelerations. The main motivation for this choice was 
the high symmetry of the problem, which allows the restriction to 
a low number of harmonics. This in turn reduces the computational 
elTort and enables a larger number A' of bodies and hence higher nu- 
merical resolution than with more traditional methods, such as the 
Barnes & Hut ( 1986) tree code. In fact, Weinberg & Katz (2002) 
even claim that cusp removal cannot be successfully simulated by 
a tree code (unless N » 10^), because the noise in the forces would 
scatter bodies off the resonant orbits, which, according to these au- 
thors, are driving the angular momentum transfer. Because of this 
argument, it was not entirely clear whether the aforementioned re- 
sults of fully self-consistent simulations are realistic or not. 

There is, however, a price to pay for using spherical harmon- 
ics. Most importantly, any method based on spherical polar coor- 
dinates is vulnerable to off-centring: if the modelled density distri- 
bution is not centred on the origin of the coordinate system, many 
terms are required for the series to converge. Conversely, when off- 
centring occurs, the gravity obtained from only the first few terms 
is strongly biased, which may induce artificial m = 1 instabilities 
and/or interfere with natural m = 1 instabilities. The only proper 
way to avoid this problem is to enforce centring, i.e. not allowing 
coefficients with odd /, m to carry any weight. Clearly, this makes 
systems with any properties that are not reflection symmetric, in- 
cluding instabilities, impossible to simulate in this way. 

IWeinberg & KatJ i2002h included coefficients of odd {l,m) 
and hence their simulations were vulnerable to artificial effects due 
to off-centring. Sellwood (2003) could reproduce their results, but 
if he suppressed the coefficients with odd (Z, m), the results were 
completely altered in that more than hundred instead of a few bar 
rotations were required to rem ove the density cusp. Sellwood also 
pointed to a previous study by White {1983) which reported that 
the I = 1 term can cause num eric al artifacts. Thus, apparently, the 
simplified bar-halo model of [Weinberg & Katz undergoes a non- 
reflection symmetric evolution. However, owing to the problems 
inherent to the use of spherical harmonics, it is unclear from the 
studies of Weinberg & Katz and Sellwood ( 2003) whether this evo- 
lution is an artifact of their algorithm (i.e. the usage of spherical 
harmonics) or not. 

In an effort to resolve this question we performed similar sim- 
ulations ourselves with a different type of A'-body code, which does 
not rely on an expansion in spherical harmonics and hence avoids 
these problems. Technical details of the models used, and the A'- 
body code are given in ^ our results are presented in ^ We dis- 
cuss our results and draw conclusions in 341 



2 MODELLING 



For the majority of our simulatio ns we model the halo as an 
isotropic spherical lHemguisJ il990l) model, which has density 



P(r)-- 



(2) 



with Miiaio and the total mass and scale radius, respectively. This 
profile has the same 'cuspy' behaviour as the NFW profile Q in 
the limit r — » 0, but has the virtue of having less mass in the outer 
regions of the halo. This means that a simulation with the same par- 
ticle number will have greater resolution in the cusp. We introduce 
the bar as an external potential applied to the A'-body simulation. 
Only the quadrupole moment of the bar is included in the force cal- 
culation, since the monopole term would have the effect of adding 
mass at the centre of the halo, altering the equilibrium, and higher 
order terms are far weaker, having little effect on the evolution, and 
are ignored for simplicity. A convenient fitting formula, which be- 
haves prop erly as a quadrupole for both r — > and r ^ oo is 
provided bv lHemguist & Weinberg! il992h : 



Oh: 



+ rl) 



sin cos 2((/) - 



(3) 



Here, r, = r/a with a the semi-major axis of the bar, r, 6, and ip are 
the usual spherical polar coordinates, and Mba, and £lp the mass and 
pattern speed of the bar, respectively, a, p and y are dimensionless 
parameters determined by best fit to the chosen bar density profile 
and shape. The bar is 'turned on' adiabatically over ~ 10 rotation 
periods to minimise transient effects. 

In th e majority of ca ses we follow iHemguist & Weinbere 
i 19921) an d ISellwooJiiooi) and model the bar as an n = 2 lFerrersl 
118771) bar with axis ratio 1 : 0.5 : 0.1. This is fitted by the pa- 
rameters cy 0. 1404, p ^ 0.4372 and y = 2 (while an incorrect 
version of this fitting formula is reproduced by Sellwood (200^), 
the correct model was used in his simulations, Sellwood, private 
communication). 

[ 1:; simulations of IWeinberg & Kat^ i2002h were of a NFW 
density profile halo 0, truncated at large radii (as the NFW model 
has infinite mass, if considered for r — » oo). The bar density 
profile used was that of a homogeneous ellipsoid of axis ratio 
1 : 0.5 : 0.05. This is well fitted by the same fitting formula 
used for the n=2 Ferrers bar with the parameters a ^ 0.1227, 
p ^ 0.6288 and y = 5.314. In Section lT2l we present simulations 
with the same parameters as i Weinberg & Katz ( 2002) in an effort 
to directly compare our simulations with theirs. For this we used a 
truncated NFW halo which we defined as having the density profile 



p(r): 



(r/r,)(l + r/rj: 



■sech(r/10rs) 



(4) 



2.1 Technical details 

We generate the initial positions and velocities from the den- 
sity j2j and the isotropic distributio n function. In t he case of 
the Hernquist halo this is known <Hernauisll Il990l) . while in 
the case of the truncate d NFW halo it can be found numer- 
ically using Eddington's Il916l formula. Quasi-random numbers 
were employed in order to suppress Poisson noise. The A'-body 
simulations were performed using the publicly available A'-body 
code gyrfalcON, which is based on Dehnen's t2000; 2002) force 
solver falcON, a tree code with mutual cell-cell interactions and 
complexity 0(N). The performance of the code is comparable 
to that of the spherical-harmonic based so-called SCF method 
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jHemauist & OstrikeJ [1993 IWeinberd [l999l) with terms up to 
n,l = S included. 

For the Hemquist halo we used units of time, mass, and length 
such that G = 1, = 11 and M^aio = 1. For the truncated NFW 
halo we used units such that G = 1, Ts = 1, and Vdicmax = 1 The 
equations of motion were integrated using the familiar leap-frog 
integrator with minimum time step 2"' and a block-step scheme 
allowing steps up to four times larger. Individual particle time steps 
were adjusted in an (almost) time- symmetric fashion such that on 
average 



T; = mm 



0.01 
1^' 



0.01 

w 



(5) 



with O, and a, the gravitational potential and acceleration of the ilh 
body. A fiducial simulation with N = 3 x l(f and no imposed bar 
conserved energy to within 0.2% over 1500 time units correspond- 
ing to ~ 120 bar rotations. 

Throughout this paper lengths will be quoted in terms of halo 
scale lengths . 



3 RESULTS 

Initial simulations of a Hernquist model halo used A' = 300 000 and 
a bar of length (semi-major axis) a = 0.7 and mass 30% of the halo 
mass interior to r = a, as it was for all simulations in this halo type. 
The bar rotates with a fixed pattern speed with co-rotation at r^,, in 
this experiment and all others. These experiments showed a rapid 
change in the 1% Lagrange radius (the radius within which 1% of 
the total mass of the halo is contained, Fig.Q. This only demon- 
strates that the density at the origin was reduced by the action of 
the bar, which is not necessarily the density peak of the halo. The 
spherically averaged density profile (Fig.|5J seems to show the for- 
mation of a core after just 8 bar rotations. Neither of these measures 
take into account any possible change in the position of the centre 
of the distribution. Similar simulations with N ranging from lO'' to 
10* showed essentially the same behaviour. This is in agreement 
with the appare nt rapid evolution seen in simulations including the 
I - I term by ISellwood i2003h . though the timescale is slig htly 
longer . It is similar to the evolution observed by Weinbere & Kata 
i2002h . though the halo model is different (see Section 13. 2> . 



3.1 Unstable evolution 

This reduction in density is, however, not demonstrative of the cre- 
ation of a cored density profile. Simply plotting the particle posi- 
tions in the x-y plane (Fig.|3j shows that within a few bar rotations 
the densest part of the halo is no longer at the origin. This creates 
the appearance of a cored density profile if analysis is performed 
with respect to the origin. This effect occurs even if the initial dis- 
tribution of bodies is made exactly symmetrical initially (for each 
body at x, v there is another one at -x, -v). 

This off-centring means that these simulation results are in- 
valid. In this simple model, the bar's centre is fixed to the origin 
and does not (is not allowed to) react to the off-centring of the halo. 
This leads to a wildly unphysical interaction between the bar and 
the halo which rapidly transports energy and angular momentum to 
the halo, and invalidates all results found this way. We strong ly sus- 
pect that the rapid evolution rep orted bytW einberg & Katj J2002h 
and its Z = 1 dependence seen bv lSellwood 42003) are actually due 
to this off-centring. Since this effect is reproduced in our simula- 
tions, it is clear that it is a real instability of the simplified bar-halo 
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Figure 1. Time evolution of the 1% Lagrange radius in flawed simulation. 
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Figure 2. Spherically averaged density profile (inner part) for the same sim- 
ulation as in Fig.0 Shown are the initial profile {dashed); and the profile 
after 8 bar rotations, considered about the original centre (solid), or about 
the centre of the cusp {dotted). The profile appears to be cored after 8 bar 
rotation times but this is because the implicit assumption of conserved (or at 
least approximately conserved) spherical symmetry is incorrect, see Fig.l3l 
The profile considered about the centre of the cusp shows that the cusp is as 
strong as ever, but the force applied by the bar has become wildly unphysi- 
cal, so this is result is, in any case, insignificant. 
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Figure 3. Positions in the x-y plane of bodies witli \z\ < 1 after 8 bar rota- 
tions. The plot shows that the centre of the distribution has shifted from the 
origin, thus invalidating the simulation setup. 



system investigated, and is not entir ely due to an art ifact of the 
spherical harmonic approach noted bv lSellwood i200.'3) . 

The initial conditions with the applied quadrupole potential 
are entirely point symmetric about the origin. However they are un- 
stable to motion of the cusp away from the centre of the bar. This 
occurs even when the initial particle distribution is defined to be ex- 
actly symmetric, because slight asymmetries caused by numerical 
truncation errors break the imposed symmetry. 

The cause of this instability can be understood by considering 
the external quadrupole potential applied to the W-body system. 
Figure shows the position of the central density in the frame of 
the rotating bar, overlaid upon a contour plot of the quadrupole 
potential. The cusp of the halo is initially at a saddle point of the 
external potential. This is unstable to a small movement of the cusp. 
We can think, in a simplistic picture, of the 'cusp' falling into a 
potential well. The restoring force due to the bulk of the halo is too 
weak to counter this. We also note that the cusp does not fall to the 
exact centre of the well, which is presumably due to the Coriolis 
force. 



3.2 NFW halo 



In an effort to re-examine the results of IWeinberg & Kat^ l2002h . 
we ran simulations of an NFW halo with the same parameters they 
used. These simulations had a short, light bar with semi-major axis 
a=0.5ri; and a mass equal to 15% of the halo mass within the bar 
radius. Since Weinberg & Katz suggested that high numerical res- 
olution was required to accurately simulate the resonant dynamics, 
we modelled the halo with particles of varying mass, with the low- 
est mass particles in the inner regions where the important resonant 
processes occur. This allows an substantial increase in the reso- 
lution of the inner regions of the NFW halo, without increasing 
the total particle number. We can define Neff, the effective parti- 
cle number, as N^^ = ^tot/™lowest' where Mjot is the total mass 
of the halo, and mjowest '■^^ mass of the lowest mass particles 
(the particles in the innermost regions of the halo). In this way we 
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Figure 4. Motion of the density centre in the plane co-rotating with the 
bar during the first 11 bar rotations (solid line) for the same simulation 
as in previous Figures. The thick dashed line indicates the orientation of 
the bar, which rotates anticlockwise. Contours of the sum of the effective 
bar quadrupole potential in this frame and an approximation for the halo- 
bulk potential (halo without cusp, using a y = Dehnen 1 1993) model) 
are shown as dotted curves. Stationary points perpendicular to the bar are 
maxima, origin is a saddle point, other stationary points along the bar axis 
are minima. 



were able to improve the resolution of our simulations such that 
^eff - 5 X N without loss of stability. 

Experiments utilising gyrfalcON, and with the N^fj- varying 
between 30 000 and 6 000 000 showed no significant evolution of 
the halo under these conditions. The bar potential was too weak 
to destabilise the halo and move its centre. This is in contrast to 
simulations with N betwee n 10 000 and 4 00 000, run with the po- 
tential expansion codes of IWeinberg & Kata. and . Sell wood which 
did become unstable rapidly (within ~10 bar rotations). The very 
low A' at which this instability is observable in the potential ex- 
pansion case, and the absence of instability even in the high A' tree 
code (gyrfalcON) case, indicates that this difference is not caused 
by particle noise drowning out the effect of the bar in the tree code 
case. The potential expansion approach is clearly more susceptible 
to the instability of the system than a tree code. 

It should be noted that there are two separate instabilities. The 
genuine instability of the artificial bar-halo system was responsible 
for the shift of the cusp away from the origin in the simulations 
with the Hernquist halo. In the NFW case the bar is too weak to 
destabilise the halo when gyrfalcON is used, but using a potential 
expansion approach exacerba tes the pro blem because of the odd 
(/, m) instability referred to by IS ell wood 



3.3 Preventing instability 

3.3.1 Forcing symmetry 

In our simplified model, both the halo and the bar are entirely point 
symmetric about the origin. The particle distribution would remain 
so for all time (thus preventing the instability) if it were not for trun- 
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cation errors in tlie properties of tlie particles as determined during 
the simulation. In order to investigate the situation without the con- 
fusion caused by the instability, we ran a number of simulations in 
which the distribution of bodies was kept point symmetric about 
the origin at all times. To this end we treat the bodies as pairs, and 
for each pair a and b, we enforce that x^, + Xj, = Q and v„ + = 
in the initial conditions and after each time step. 

This is effectively the equivalent of simulations using an ex- 
pansion in spherical harmonics which exclude all the terms wit h 
odd (/, m). Such experiments were carried out bv lSellwoodj2003h . 

As Figure |5| shows, the ability of a rotating bar to remove 
a cusp is genuine, and not reliant upon instability. Even when 
the simulation is completely symmetric, angular momentum trans- 
ferred from the bar to the halo eventually removes the cusp. The 
dependence of the effect on A' and bar strength are comparable to 
the results obtained by Sellwood (Figure|5j. The major difference 
to note is that our results do not yet seem to converge at high A'. 
This is likely caused by our force solver allowing stochastic jitter 
in the body distribution to be passed to the forces more readily than 
a method involving only low-order spherical harmonics does. 



3.3.2 Using a more realistic bar model 

An alternative method for stabilising the simulation is to allow the 
bar to "have mass", and thus provide a restoring force on the halo. 
This was achieved in our simulations by including the monopole 
component of the bar potential in the force calculation. The ini- 
tial distribution function of the halo with the additional monopole 
potential can be determined using Eddington's formula. In these 
simulations the monopole component of the potential was included 
at all times, the quadrupole was added adiabatically as in all pre- 
vious simulations. Any off-centring of the halo causes an artificial 
transfer of linear momentum from the bar to the halo, so in addi- 
tion the bar was constrained to move in such a way that the linear 
momentum of the bar-halo system was conserved. 

Figure |5| shows a comparison of the results from this ap- 
proach to those from the symmetrised simulations. The movement 
required of the bar in order to conserve the linear momentum of 
the system was extremely small, and only carried the bar a dis- 
tance ~10"^r,, which is several orders of magnitude smaller than 
the movement of the cusp seen in Figure|4] Simulations which in- 
cluded the monopole in the bar potential, but in which the bar was 
piimed to the origin produced results which were nearly indistin- 
guishable from those with a moving bar. It is clear that the restor- 
ing force provided by inclusion of the monopole term is sufficient 
to prevent the growth of the instability in this case. The presence 
of the monopole and removal of the artificially enforced symmetry 
prevent the rapid "runaway" evolution seen at f ^ 170 in the cor- 
responding symmetrised simulation. However the bar still removes 
the cusp of the halo over approximately the same timescale in the 
two simulations. 



3.4 Angular Momentum conservation 

These experiments are deliberately kept as simple as possible, 
which leaves them unrealistic. One major difference between the 
simulations described thus far and the true situation is the apparent 
infinite supply of angular momentum for the bar, which keeps it 
rotating with the same pattern speed for the entirety of the simula- 
tion, however much angular m omentum has been transferred to the 
halo. IWeinberg & KatJ i2002l) pointed to a "suite of simulations 
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Figure 5. The 1 % Lagrange radius plotted against time for many different 
simulations. Top: the effect of increasing bar size with fixed N = 300 000 
and enforced symmetry. Bar mass was defined as 30% of the halo mass 
within the bar radius for all cases. Also shown are the results of a fiducial 
run with N = 300000 and no imposed bar. Bottom: the effect of increasing 
resolution (paiticle number), for fixed bar size a = 0.7 and enforced sym- 
metry also showing (Dotted) simulation with conserved lineal' momentum 
with same bar size, and A? = 1 000 000 (*) 



with a slowing bar, whose pattern speed follows from the conser- 
vation of angular momentum of the combined bar-halo system" as 
further evidence pointing to the effectiveness of this process in re- 
moving the cusp. However these results are tainted by the instabil- 
ity of the system, so we repeat the experiments with a slowing bar 
and the symmetrisation method to prevent instability. This is some- 
what similar to experiments performed by Sellwood 1 2003), with 
the / = 1 spherical harmonic suppressed. Since initially the bar has 
zero mass, and the mass (and thus angular momentum) is increased 
adiabatically, angular momentum is introduced to the system by the 
growing bar. We calculate the angular momentum of the halo after 
each time step and use that to calculate the change in the pattern 
speed of the bar at every point, so that the total angular momentum 
added to the halo is no more than that of the bar at its full mass, 
rotating with its original pattern speed. 

As Figure|6|shows, when the bar is not given an infinite supply 
of angular momentum, the situation is somewhat different. Essen- 
tially all the angular momentum of the bar is absorbed by the halo 
within 15 initial rotation periods (in fact the bar has made ten com- 
plete revolutions at that point). This absorbed angular momentum 
is insufficient to remove matter from the cusp to any great extent. 
In experiments carried out without symmetrisation, the central den- 
sity did fall to approximately that seen in the experiments with a bar 
rotating with a constant pattern speed. This was again because the 
cusp of the distribution shifted from the origin. 



4 CONCLUSIONS 

We find that over-simplifying the interaction of a barred galaxy 
with a cusped CDM halo may lead to substantial over-estimation 
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Figure 6. Density profile for a slowing bar m a N = 300 000 representa- 
tion of a Hemquist halo initially {dashed) and after 16 initial bar rotation 
periods (solid). The bar has not gone through 16 rotations, because it has 
decelerated, and by this point is rotating at 10% of its initial pattern speed. 



of the reduction in central density. The cause is an instability of 
the simplified model: a shght off-centring of the CDM cusp from 
the origin results in a net force from the bar, which is fixed to stay 
centred, further driving the cusp away from the origin. The insta- 
bility is implicit in the model and is not simply a result of deter- 
mining the potential as a spherical harmonic expansion about the 
origin. It seems very lik ely that the rapid evolution reported by 
IWeinberg & Katd j2002h and Sellwood 1 2003, when odd / terms 
where included in the computation of the forces) were an artifact of 
this instability in conjunction with averaging over spherical shells 
concentric with the origin (rather than the density centre). 

If this artificial instability is suppressed, either by enforcing 
reflection symmetry w.r.t. the origin, or by including the monopole 
of the bar potential and conserving linear momentum, the removal 
of the CDM cusp is a very slow process, requiring (9(100) bar- 
rotations, i.e. a significant fraction of a Hubble time. Moreover, 
this simplistic model contains an additional artifact: the infinite 
supply of angular momentum from the ever torquing bar. If this 
is removed by controlling the total angular momentum, the cusp is 
at most slightly modified. This means that, under these conditions, 
the amount of angular momentum needed to remove a CDM cusp 
is far larger than that present in the inner parts of disk galaxies. 

Of course, this work makes a large number of assumptions. 
The halo has been assumed to be spherically symmetric, isotropic, 
without substructure, isolated and initially in equilibrium. The bar 
was assumed to be rigid, and gas physics and bar-disk interaction 
were neglected. These are the same assumptions made by the paper 
that proposed this mechanism ( Weinberg & Katz 2002). These sim- 
ulations cannot tell the whole story, but they do provide an insight 
into the fundamental physical mechanism. 

We conclude that, contrary to the original proposal, these 
simulations suggest that angular momentum (and energy) trans- 
port from a disk-bar is not an effective way to destroy CDM 
cusps, in agreement with studies using less simplistic models 



JValenzuela & KlvpinlEo03LlAthanassoulall2004 . If galactic scale 
CDM halos truly are cored, it seems that some other process is re- 
sponsible. 



ACKNOWLEDGEMENT 

The authors would like to thank their referee, Kelly HoUey- 
Bockelmann, for her constructive suggestions. PJM acknowledges 
the support of the UK Particle Physics and Astronomy Research 
Council (PPARC) through a research student fellowship. Astro- 
physics research at the University of Leicester is also supported 
through a PPARC rolling grant. 



REFERENCES 

Athanassoula E., 2004, in lAU Symposium Vol. 220, Bars and the 

connection between dark and visible matter, p. 255 
Barnes J., Hut P, 1986, Nature, 324, 446 
Binney J., Gerhard O., Silk I., 2001, MNRAS, 321, 471 
Binney J. J., Evans N. W., 2001, MNRAS, 327, L27 
Blais-Ouellette S., Amram P, Carignan C, 2001, AJ, 121, 1952 
Blumenthal G. R., Faber S. M., Flores R., Primack J. R., 1986, 
ApJ, 301, 27 

Bode P, Ostriker I. P, Turok N., 2001, ApJ, 556, 93 
Cote S., Carignan C, Freeman K. C, 2000, AJ, 120, 3027 
Coh'n P, Avila-Reese V., Valenzuela O., 2000, ApJ, 542, 622 
de Blok, W. J. G., McGaugh, S. S., & Rubin, V. C. 2001a, AJ, 122, 
2396 

de Blok, W. J. G., McGaugh, S. S., Bosma, A., & Rubin, V. C. 

2001b, ApJl, 552, L23 
Dehnen W., 1993, MNRAS, 265, 250 
Dehnen W., 2000, ApJl, 536, L39 
Dehnen W., 2002, J. Comp. Phys., 179, 27 
Eddington A. S., 1916, MNRAS, 76, 572 
El-Zant A., Shlosman I., Hoff^man Y., 2001, ApJ, 560, 636 
Ferrers N. M., 1877, Q. J. Pure Appl. Math., 14, 1 
Gentile G., Salucci P, Klein U., Vergani D., Kalberia P, 2004, 

MNRAS, 351, 903 
Hayashi E., Navarro J. F, Power C, Jenkins A. R., Frenk C. S., 

White S. D. M., Springel V., Stadel J., Quinn T. R., 2004, MN- 
RAS, 355, 794 
Hernquist L., 1990, ApJ, 356, 359 
Hernquist L., Ostriker J. P, 1992, ApJ, 386, 375 
Hernquist L., Weinberg M. D., 1992, ApJ, 400, 80 
Holley-Bockelmann K., Weinberg M. D., Katz N., 2003 

astro-ph/0306374 
Kaplinghat M., Knox L., Turner M. S., 2000, Physical Review 

Letters, 85, 3335 
Moore B., Governato F, Quinn T, Stadel J., Lake G., 1998, ApJl, 

499, L5-^ 

Navarro J. F, Frenk C. S., White S. D. M., 1996, ApJ, 462, 563 
Power C, Navarro J. F, Jenkins A., Frenk C. S., White S. D. M., 

Springel V., Stadel J., Quinn T., 2003, MNRAS, 338, 14 
Salucci P, Walter F, Borriello A., 2003, Aap, 409, 53 
Sellwood J. A., 2003, ApJ, 587, 638 
Sellwood J. A., 2005, ApJ, submitted <astro-ph/04075^ 
Spergel D. N., Steinhardt P. J., 2000, Physical Review Letters, 84, 

3760 

Swaters R. A., Madore B. F, van den Bosch F. C, Balcells M., 
2003, ApJ, 583, 732 



Halo evolution in the presence of a disc bar 7 

Trott C. M., Webster R. L., 2002, MNRAS, 334, 621 

Valenzuela O., Klypin A., 2003, MNRAS, 345, 406 

van den Bosch F. C, Swaters R. A., 2001, MNRAS, 325, 1017 

Weinberg M. D., 2004, astro-ph/0404169 

Weinberg M. D., 1999, AJ, 117, 629 

Weinberg M. D., Katz N., 2002, ApJ, 580, 627 

White S. D. M., 1983, ApJ, 274, 53 



